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The metallic surface states of a topological insulator support helical Dirac fermions protected by 
topology with their spin locked perpendicular to their momentum. They can acquire mass through 
magnetic doping or through hybridization of states on opposite faces of a thin sample. In this case 
there can be a component of electron spin oriented perpendicular to the surface plane. The electron- 
phonon interaction renormalizes the dynamics of the charge carriers through their spectral density. 

It also modifies the gap channel and a second spectral function enters which, not only determines 
the out of plane spin component, but also comes into in-plane properties. While the out of plane 
spin component is decreased below the Fermi momentum (fcf), the in plane component increases. 

There are also correlation tails extending well beyond k,F- The angular resolved photo-emission 
line shapes aquire Holstein side bands. The effective gap in the density of states is reduced and 
the optical conductivity aquires distinct measurable phonon structure even for modest value of the 
electron-phonon coupling. 

PACS numbers: 75.70 Tj, 71.38 Cn, 78.67.-n 


I. INTRODUCTION 

Helical Dirac fermions protected by topology and 
found at the surface of topological insulators (TI)i~4 ex¬ 
hibit spin locked perpendicular to their momentum.^ 8 
Doping with magnetic impurities can break time rever¬ 
sal symmetry and create massive Dirac fermions as has 
been experimentally demonstrated^ in B^Se^. A gap 
can also be introduced in thin films when the distance 
between top and bottom surfaces is of the order of the 
extend in space of the surface statesJ4iL Gapped Dirac 
fermions arise in many other systems, for example in 
two dimensional membranes such as the dichalcoginide 
MoS- 2 ^— and silicenei^r— with buckled honeycomb lat¬ 
tice. In both these cases pseudospin plays the analogous 
role to the real spin of topological insulators. 

The electron-phonon interaction renormalizes quasi¬ 
particle dynamics and leads to important observable 
changes in electronic properties 22-3 ^ which illustrate the 
effects of many body renormalizations not captured in 
single particle theories. For the Dirac electrons in 
graphene as an example, features observed in the den¬ 
sity of states^ - — and in the dispersion curves measured 
in angular resolved photo emission spectroscopy^ have 
been interpreted as phonon structure. 36 The optical prop¬ 
erties of graphene are also renormalized in a non trivial 
way j 22 i 25 i 26 i 37 ' 38 In a simple bare band picture there is no 
optical absorption in the photon energy region between 
the Drude intraband contribution centered about w = 0 
and the interband onset at twice the value of the chemical 
potential w = 2[i. In reality the real part of the dynamic 
longitudinal conductivity is observed to be finite and al¬ 
most one third its universal backgroun d 39 ' 40 value in this 
photon region. 3 ' This absorption is due to many body 
renormalizations and is at least partially assigned^ 2 ’™ 


to the electron-phonon interaction which provides boson 
assisted processes referred to as Holstein processes. A 
phonon is created along with an electron-hole pair. 

Understanding the transport properties of the Dirac 
electrons on the surface of a topological insulator is im¬ 
portant for possible device applications. At finite tem¬ 
perature the electron-phonon interaction is expected to 
be an important scattering chanel^L, X. Zhu et al42 
have studied the surface phonons on the (001) surface 
of B^Sej, and in particular have found a giant Kohn 
anomaly associated with a branch having a maximum of 
18 THz. From measurements of the phonon self energy, 
the same group^ determined the size of the electron- 
phonon interaction report a coupling constant of 0.43 for 
a particular branch which is much larger than reported in 
angular resolved photo-emission spectroscopy. One such 
study by Z. H. Pan et al4i gives a mass enhancement A of 
0.08 while another by R. L. Hatch et al4I found A ~ 0.25. 
A more recent report^ gives ~ 0.17 with characteristic 
phonon energy ~18 meV. 

The results of X. Zhu et al. are further supported by 
an infrared study by A. D. LaForge et al4£ which found 
a strong electron-phonon coupling to a 7.6 meV optical 
phonon while S. Giraud et al. 47 provide arguments for 
coupling to acoustic phonon with A ~ 0.42 in their films 
and even larger in other geometries 4^ Very recently J. A. 
Sobota et alM, using tinre-resolved photoemission spec¬ 
troscopy, find evidence for coupling to a 2.05 ThZ surface 
phonon mode with the Dirac electrons in B^Se^. Possi¬ 
ble complications in the interpretation of optical pump¬ 
ing on time resolved data were discussed by S. Ulstrup^I 
and need to be kept in mind. Finally, the recent transport 
measurements of M.V. Costache et al4L were interpreted 
with strong coupling to a single optical phonon mode of 
energy ~ 8 meV. 

In this paper we will emphasize the effect of the 
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FIG. 1. (Color online) Illustrative representation of Dirac 
cones (left frame) and including a gap A (right frame). 



electron-phonon interaction on the spin texture of gapped 
Dirac helical fermions. We will present for comparison, 
results for the familiar spectral density {A I (k,u)) associ¬ 
ated with quasiparticle renormalization which we com¬ 
pared with the much less familiar equivalent function 
which enters gap renormalizations ( A z (k , w)). This is the 
essential quantity for the calculation of the z-axis (per¬ 
pendicular to the surface states) properties. While only 
Ai(k , u) is needed in calculations of the density of states 
N(lo) both A z and A/ enter the dynamic conductivity 
Vxx(u). In section II we present the necessary formal¬ 
ism for both the self energy associated with quasiparti¬ 
cle renormalizations £/(u;)) and with the renormalization 
of the gap £ 2 (w)) which determine the needed spectral 
functions A/(fc,w) and A z (k,u>). Numerical results are 
presented in section III which deals with spectral densi¬ 
ties, density of states and optical conductivity. In section 
IV we consider both in plane and out of plane spin tex¬ 
ture. A summary and conclusions are given in section 
V. 


II. FORMALISM 

We begin with a model hamiltonian for the electronic 
states at the surface of a topological insulator which has 
the form 


H 0 = hv F [k x a y - k y cr x \ + Aa z + I?o(k) (1) 

where cr’s are Pauli spin matrices for real electron spin, 
v F is the Fermi velocity of the Dirac electrons, k is mo¬ 
mentum, A is a gap which can be introduced for exam¬ 
ple by doping with magnetic impurities and I?o(k) is a 
quadratic term A’o(k) = with m the electron mass. 
The introduction of Eq gives particle hole asymmetry. 


FIG. 2. (Color online) The real and imaginary part of the 
electron-phonon self energies E J (ta) and E z (ca) in meV as 
a function of t j in meV. E z (ca) (real part, solid black), and 
imaginary part (dashed red) gives the quasiparticle renormal¬ 
ization and E z (a>) (real part, dotted blue), and imaginary 
part (dash-dotted purple) gives the gap renormalization. Our 
choice of g=10 corresponds to a mass enhancement A ~ 0.3. 


Here we will assume that m is very large, as a first ap¬ 
proximation, and will ignore this term. 

The non interacting Green’s function takes on the form 


G 0 (k, s,iu n ) 


1 

ico n + n — Sy/h 2 Vpk 2 + A 2 


( 2 ) 


with iu> n the Matsubara frequencies and /r is the non 
interacting chemical potential. The eigen energies are 
£k, s = s^/h 2 Vpk 2 + A 2 with s = ± for conduction and 
valence band respectively. These dispersion curves are 
shown schematically in Fig. 1 with (right frame) and 
without (left frame) a gap. 

The matrix (2 x 2) Green’s function Go(k, iuj n ) can be 
written in terms of the scalar Go(k, s, iui n ) as 


G 0 (k,iui n ) 


\ ^(1 + sFk ■ o-)G 0 (k, s, iuj n ) (3) 


with the vector F& defined as 

_ {-hv F k y ,hv F k x , A) ^ 

k ~~ y/H 2 v 2 pk 2 + A 2 

We want to include an electron-phonon interaction. The 
simplest is the Holstein model for coupling to a phonon 
mode of energy lo f with matrix element between elec¬ 
tronic and phonon assumed to be a constant (g). It is 
written as 

H e —ph — g^E ^ ^ c k,s^k',s(bj i /_k V k 7 ) (5) 

k.k',s 
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where c^ g creates an electron of momentum k and spin 

s and &£,_ k creates a phonon of momentum k' k re¬ 
duced to the first surface state Brillouin zone and energy 
loe with g coupling constant. The perturbing Hamilto¬ 
nian (5) provides a self energy to electron motion which 
has two components one proportional to the unit matrix 
I and the other to a z . Writing E(iw n ) = E J (iu n )I + 
E z (im n )a z we get in lowest order perturbation theory 


E 1 (iu> r 


) = 


g 2 U> 2 E 

2 ^ 

q,s 


/f(%«) + /b(wb) . /b(wb) + 1 - /j?(e ? , s ) 


+ H + WE - £q, s iu n + /J ~ U) E ~ £, 


q,s 


( 6 ) 


and 


E z {iu> n ) = 


2 2 

= ^ 
q.s 


sA 


2 \Zh 2 Vpq 2 + A 2 

fF{s q ,s) + Ib{^e) Ib{^e) + 1 ~ f F {Sq,s) 


+ p, + wg - + n - u} E ~ £, 


q,s 


( 7 ) 


where f F and fs are fermion and boson distribution 
functions l/[e^ u ~^^ kBT ± 1] respectively. Here g is the 
chemical potential which applies only to the electron dis¬ 
tribution. Note that E z (ioj n ) in Eq. (7) is directly pro¬ 
portional to the gap A and will vanish for A = 0. In 
terms of the self energies T, 1 (ioj n ) and E z (tw n ) the inter¬ 
acting matrix Green’s function takes on the form 


G(k, iuj n ) = i '^2 (1 + sH k • a)G(k , s, iu n ) 

s=± 


with 


H 


k = 


( - hv F k y , hv F k x , A + E z {im n )) 
K 2 v 2 F k 2 + (A + E z (im n )) 2 


( 8 ) 


( 9 ) 


and 

G(k, s, iui n ) = 


iu n + fi — E 7 (iw n ) — Sv7(A"+ E z (iuj n )) 2 + h 2 v 2 F k ^ ^ 

In Ref. (19) a factor E z *(zw n ) was mistakenly introduced 
instead of E z (iw n ) in (9) and (10). This leads to small 
numerical differences but has no qualitative significance. 
The spectral function associated with / and a z matrix 
are 

Aj(k , s, ui) = - ImG(k , s, iu n —> oo + id) 

7T 

= A{k,s,m) (11) 

which is the quantity measured in angular resolved photo 
emission spectroscopy ARPES^ Further 

A z (k,s,m) 

= - ImTr[a z G{k , s, iu n —> to + z<5)] 


1 , s[A + E z (w + id)]G(k, s, uj + id) 


= - Im{ 


f \Jh 2 v‘pk 2 + (A + E z (lu + id)) 2 


}. ( 12 ) 


A considerable and instructive mathematical simplifica¬ 
tion of these complicated expressions for (11) and (12) 
result when the imaginary part of the z-component of 
the self energy E z (*w n —> u + id) is ignored. We get 


A(k, s,u>) = — 


ImT, T (u)) 


(13) 


with 


and 


f [oj — sVM] 2 + [ImE f (w)] 2 
A z (k,s,uj ) 

s[A + ReT, z (oj)]A(k, s, m) 

\J H 2 Vpk 2 + (A + Re E z (w)) 2 

u; = u> + /j, — ReE 1 {u>) (15) 


(14) 


M = h 2 v 2 pk 2 + (A + ReE z {co)y 


(16) 


where M is a function of k and to. If we further take 
the imaginary part of E / (w) to be zero (13) reduces to 
a Dirac delta function but the real part of both E J and 
E z remain which renormalize the single particle energies 
and the gap. The probability of occupation of the state 
k denoted by n(k) at temperature T is given by22 

/ +oo 

dwA(k, s, ui)f(u>) (17) 

-OO 

and this is to be compared with the corresponding ex¬ 
pression for the z-component of spin S z ( k) 


/ +oo 

du>A z (k,s,u)f(cj). 

-OO 


(18) 


We can also calculate the average value of the square root 
of the sum of the squares of x and y component of spin 
which remains locked perpendicular to momentum but its 
magnitude is changed by the electron-phonon coupling 

Js 2 (k) + S 2 (k) = I dujA x _ y (k, s, m)f (w) (19) 

J —OO 


with 

A 1 r [hvFkjGjk^^) n 

A x - y (k, s,ui) 

( 20 ) 

The density of electronic states N(u>) follows from 
A(k, s, w) on integration over k 


N (u) = Y2 A (k,s,uj). 


( 21 ) 


This quantity enters scanning tunneling microscopy ex¬ 
periments (STM) i22r— The real part of the dynamic lon¬ 
gitudinal conductivity which gives the absorption spec¬ 
trum for light follows as^ 


Rea xx (u) = 


a 2„,2 ,_2 


■E 


dm 


m ^.l-oa 2f 


[/(w) - /(w + H)] 


x[A/(k, w)A/(k,w + Q) — A z (k, uj)A z {k,m + fl)] (22) 
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with Aj( k,w) = A/(fc, s, w) and A z ( k, w) = 
A z (k, s,lo). Note that both A z and Aj enter this 
quantity even though it is an in-plane property. 


III. NUMERICAL RESULTS FOR SELF 
ENERGY AND DERIVED QUANTITIES 

In lowest order perturbation theory the chemical po¬ 
tential g in Eq. (6) and (7) is to be interpreted as hav¬ 
ing its non interacting, independent particle value and 
will be quoted in all the results to be presented in this 
paper. With electron-phonon coupling the chemical po¬ 
tential will in general be shifted from its free band value 
to Hinteracting = A* + ReT, 1 (lLO n ->• w + iS). In Fig. 2 we 
show our results for £ [iu> n —> u+iS) in meV as a function 
of energy uj in meV where we have analytically contin¬ 
ued from Matsubara to real frequencies. The solid black 
curve is Re E 4 with JmE 4 given by the dashed red curve 
while the dotted blue is for i?eE z (w) with its imaginary 
part /m£ z (u) given by the purple dash-dotted curve. 
First the real part of the quasiparticle self energy is not 
zero at oj ~ 0 but rather is of the order of ~ —5meV 
which reduces the chemical potential of the interacting 
picture by about 1 /4 of its bare band value here taken to 
be 25 meV with the gap A equal to 20 meV. With g set 
to a value of 10 in Eq. ( 6 ) and (7), the mass enhancement 
parameter A ~ 0.3 as can be seen in Fig. 2. The param¬ 
eter A is the slope of the real part of E 1 {oj) at oj = 0, i.e. 
A = {d/doj)ReT, 1 {oj) at oj = 0. We have also taken the 
Einstein phonon energy to be 7.5 meV which is represen¬ 
tative of what is seen in many but not all experiments 
in B^Se 3 and has been identified in the work of X. Zhu 
et al— as a surface phonon. As discussed in the intro¬ 
duction coupling to specific optical phonons have been 
found in other experiments^^—. Of course there can 
also be coupling to acoustic phonons as studied by Giraud 
et al.i 47 i 48 These authors use a continuum model for the 
phonons and a deformation potential model to describe 
their coupling to the helical surface Dirac fermions with 
due attention to the modifications brought about by he- 
licity of the charge carriers. In our simplified model these 
complications enter only in determining the size of g in 
Eq. ( 6 ) and (7) which we set through consideration of 
experiment (i.e. A ~ 0.3). In principle these could also 
be coupling to electronic modes but this is not considered 
in this work which deals with electron-phonon coupling 
alone. Our choice of gap A = 10 mev is representative 
of the work of Chen et al£ who found in ARPES ex¬ 
periments a gap of — — 7 mev in B^Se^ magnetically 
dopped with Mn. We note prominent phonon structure 
at w = ±wg = ±7.5 meV by choice, with additional 
structures at —g+A — oje and at — g— A — oje- For oj > 0 
the imaginary part of the quasiparticle self energy is zero 
in the interval ( 0 , teg) after which its absolute magnitude 
is finite and increases linearly with increasing value of oj. 
This dependence reflects the linear dependence on energy 
of the underlying bare density of states. Also it needs to 




FIG. 3. (Color online) The spectral density A{k,w ) as a 
function of energy oj in meV for three values of momentum 
k namely k = 0 (solid black), k = 0.01 (dotted blue) and 
k = 0.02 (dashed red) in units of the inverse lattice param¬ 
eter a. The top frame is for A z {k, oj) (associated with the 
gap renormalization) and the bottom frame is for Ai{k,oj) 
associated with quasiparticle renormalization. Vertical dot¬ 
ted black lines are at u> = ±oje with phonon Einstein energy 
set at 7.5 meV. The gap A = 20 meV and the bare chemical 
potential is g = 25 meV. 


be negative because —ImE 1 {ioj n —> ui+iS) is a scattering 
rate, and hence positive. By contrast the imaginary part 
of th e gap self ene rgy is flat because the extra factor of 
A/\A 2 u|V + A 2 appearing in (7) but not in ( 6 ) com¬ 
pensates for the density of states variarion which comes 
from the sum over q. Note also that both self energies 
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Az, T=0.1K 
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FIG. 4. (Color online) A color plot of the gap spectral den¬ 
sity A z (k,co) as a function of energy to in (eV) vertical axis 
and momentum k in (1/a) horizontal axis, with a the lattice 
parameter. The gap A = 10 meV and the Einstein phonon 
energy loe = 7.5 meV. 



FIG. 6. (Color online) The real part of the longitudinal dy¬ 
namic optical conductivity Rea xx (to) in units of e 2 /fi as a 
function of photon energy to in meV. The dotted red line in¬ 
cludes the electron-phonon interaction with g = 10 in Eq. 
(6)-(7). The solid black line, which is for comparison, is for 
bare bands (g = 0). The top frame has no gap A = 0 and the 
bottom frame includes a gap A = 20 meV. In both cases the 
bare chemical potential g = 35 meV. The arrow identifies the 
phonon structure. 



FIG. 5. (Color online) The density of states N(lo) in inverse 
eV versus energy to in eV. The dotted red line includes the 
electron-phonon interaction with g = 10 in Eq. (6)-(7). The 
solid black line which is for comparison is for bare bands (g = 
0). The top frame has a gap A = 10 meV while the bottom 
frame has a gap of 20 meV. 


have a finite imaginary part in the interval ui € (0, —uje) 
and to < — g — A — u>e- While the real part of the quasi¬ 
particle self energy can be both positive or negative, the 
real part of the gap self energy is everywhere negative 
and so decreases the bare gap A at all u considered here. 
Its magnitude is everywhere of order 1/4 of the input 
value of A with some to dependence and is encoded with 
boson structure. These self energies have a profound ef¬ 
fect on the corresponding spectral densities of Eq. (11) 
and (12) as we see in Fig. 3 where the top frame applies 
to A z (k,uj) = A z (k,s > 0,w) (conduction band) and the 
bottom frame is for Aj(k,ui) = Ai(k,s > 0, w). For the 
range of values of k shown namely, k = 0.0, 0.01 and 0.02 
in units of the inverse of the lattice parameter (a), these 
two functions do not differ much from each other. This 
is traced to the fact that A z (k,ui) has the extra factor 
, as compared with Ai(k,ui) where we have, 

,/ft 2 Upfc 2 +A 2 

for simplicity, neglected renormalizations as these do not 
change the argument in an important way. As long as 
hvFk is small compared with A the relevant additional 
factor is near one. For large k however, larger than those 
shown in Fig. 3, A z (k,w) will become much smaller in 
magnitude than Aj(k,co) by a factor of A/(hvpk). Re¬ 
turning to the top frame of Fig. 3 we first note in the 
solid black curve for k = 0 (which is the momentum 
of the bottom of the bare conduction band) a large, only 
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slightly broaden and shifted in energy, quasiparticle peak 
centered approximately at ~ 6 meV below the Fermi sur¬ 
face at lj ~ 0. For the bare bands it would be at 5 meV 
instead. But with correlations there are also additional 
features. There is a large increase in A z (k,oj) starting at 
u) = —oje = — 7.5meV followed by a large boson struc¬ 
ture extending down to — —12.5meV where it drops to 
zero. Such a feature is also seen in the dashed red curve 
but it is much smaller. The lower cut off is identified to 
correspond to the energy — /z — u>e + A. If we had shown 
the valence band contribution to the spectral density it 
would have a further boson structure at — /z — u>e — A as 
in the self energy of Fig. 2. It is clear from this descrip¬ 
tion that the effect of the electron-phonon interaction on 
the bare band is much more complicated than a simple 
constant shift in gap value and a slight shift in quasipar¬ 
ticle energies with small broadening. Note also the clear 
phonon sidebands at energies above to = to e in all three 
curves. The black vertical dotted lines identify energies 
ui = ±ug for easy reference. A different representation 
of these changes is given in Fig. 4 where we show a false 
color plot of A z {k , ui) as a function for energy ui along the 
vertical axis and momentum k along the horizontal. As 
seen in Eq. (12) and (14), A z carries the sign of s and is 
positive (red) in the renormalized conduction band and 
negative (blue) in the renormalized valence band. The 
phonon at ui = ±ub is clearly identified as is the bottom 
of the renormalized conduction band at —fi — oje + A 
and the valence band phonon structure at —/z — cue — A. 
Particularly striking is the large modification of the bare 
band dispersion curves in the region just above the renor¬ 
malized conduction band minimum. 

The large renormalization effects seen in Fig. 3 and 
Fig. 4 can have a profound effect on certain quantities 
while at the same time have much more modest manifes¬ 
tations in others as we will now describe. In Fig. 5 we 
show results for the density of states (DOS) N(u) given 
in Eq. (21). The top frame is for a gap A = 10 meV 
and the bottom for A = 20 meV. The red dotted curve 
is the renormalized quasiparticle DOS while the black 
solid curve is the bare band case shown for comparison. 
In both frames the bare chemical potential is set at 25 
meV, which is 5 meV above the gap in the conduction 
band (for the case A = 20 meV). The energy variable 
w is set so that oj = 0 corresponds to the Fermi sur¬ 
face so that the bottom of the conduction band is at —5 
meV in bottom frame and —15 meV in top frame with 
the top of the valence band at —45 meV and —35 meV 
respectively. In both cases the remaining gap between va¬ 
lence and conduction band has been very much reduced 
over its bare band value. The four phonon structures at 
ui = iwg, — /z — uje — A and — /z — oje + A are clearly seen 
with the bottom of the conduction band given by the en¬ 
ergy —fi — oje + A. By comparison, the corresponding 
boson structures in the real part of the dynamic longitu¬ 
dinal optical conductivity Rea xx (u>) which is more closely 
related to a convolution of two DOS factors and is given 
in Eq. (22), are much more modest as seen in Fig. 6. 


The top frame is for gapless Dirac fermions (A = 0) 
and is included for comparison while the bottom frame 
is for the gapped case with A = 20 meV. The solid black 
curves are for the bare band case with a small residual 
scattering rate included to broaden out the Drude peak 
due to intraband transitions. This peak is large only at 
small ui and is centered at u> = 0. The bare chemical 
potential is 35 meV and we see the onset of a second 
absorption band coming from the interband transitions 
which start at to = 2/z. These transitions continue up to 
large energies and provide the so called universal back¬ 
ground. In our units for Rea xx {uj) which is e 2 /7i, this 
backgroun d^— has a height of 1/16. As is clear in the 
top frame for A = 0 this height is almost unaffected by 
the electron-phonon interaction (red dotted curve). Stay¬ 
ing with the top frame we see however the appearance of 
the Holstein processes above w = u>e which provides sig¬ 
nificant phonon assisted absorption in the photon region 
above loe and below the main interband absorption edge 
at 2/z. The other feature to be noted is that, for the cor¬ 
related case, the onset of the interband transitions has 
moved to lower energies and is now at twice the value 
of the interacting chemical potential. The lower frame 
which is for the gap fermion case has another important 
element. As is well known ) 19 ' 39 ' 40 when A/0 there is a 
peak in the interband transitions just above the threshold 
energy which persists up to a few A above the threshold 
before the value of the background is reestablished at its 
universal value. In the clean limit we have the analytic 
result 

Rea xx (u,) = - A ") + 

/ .2 _I_ A A 2 

Ip 0 ( u — 2max(|/z|, A))] (23) 

where 9 is a heaviside function. We see that, just above 

the interband onset the conductivity is larger than its 
2 

universal value of as is seen most clearly in the solid 
black curve of the lower frame of Fig. 6 which is the bare 
band result. The horizontal straight line segment indi¬ 
cate 1/16. We note that, with electron-phonon (dotted 
red curve), the magnitude of the absorption in the region 
of the edge still remains above the universal background 
value but now there is also a small phonon structure high¬ 
lighted by the vertical arrow. No such structure is seen in 
the top frame for A = 0. It is the variation with energy 
of the background (in the presence of a finite gap) which 
allows for the phonon structure to be revealed. 

IV. SPIN TEXTURE 

The z-axis spin texture is determined from the spectral 
density A z (k,s,oj) of Eq. (12). Results are presented in 
Fig. 7 for S z (k) as a function of momentum k normalized 
to its value at the Fermi surface (fc p). S~(k) involves an 
integral over frequency ui of the overlap of A z {k , s, ui) and 
the Fermi function f(ui) as given in Eq. (18). We only 
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FIG. 7. (Color online)The z-component of spin S z (k) as a 
function of momentum k normalized to &f. The solid black 
line is for g = 0 (no electron-phonon interaction) and is for 
comparison with the dotted blue curve which includes the self 
energy with coupling g = 10 (Eq. (6)-(7)). 


FIG. 9. (Color online) Comparison of the z-component of spin 
S z (k) of Eq. (18) as a function of momentum k normalized 
to its value at the Fermi surface Icf (short dashed red curve) 
with the probability of occupation of state n(k ) of Eq. (17) 
(long dashed black curve). The heavy solid black line gives 
the case of bare bands for comparison. 


-g=0, A=20 meV 



k/k F 

FIG. 8. (Color online) The magnitude of the in plane spin 
Sx-y(k) as a function of momentum k normalized to its value 
at the Fermi surface /cf. The solid black curve is for compar¬ 
ison and involves the bare bands (g = 0) with gap A = 20 
meV. The dotted blue has the same gap but includes the 
electron-phonon renormalizations with coupling g = 10 in Eq. 
(6)-(7). The dashed red is for a smaller value of gap A = 10 
meV. 


show results for the conduction band s = +. For the 
valence band s = — 1 and S z (k) will change sign. The 
gap A = 20 meV and the chemical potential /.i = 25 meV 
as before in section III. The solid black curve is the bare 
band result for temperature T = 0.1 K. We see an abrupt 
drop to zero at k/kF = 1 because of the thermal factor 
which here is basically a Heaviside function with cut off 
at the Fermi surface. By contrast the dotted blue curve 
includes the electron-phonon renormalization discussed 
at length in the previous section. Now the magnitude 
of S z (k) at k = 0 is reduced below its bare band value. 
More importantly the jump at k/kF is smaller than it is 
in the solid curve and further, there are finite tails beyond 
this momentum which are entirely due to many body cor¬ 
relation effects that go beyond a bare band description. 
They can be taken as representative of other correlation 
effect such as those due to electron-electron interactions 
rather than electron-phonon. Compared to the modifica¬ 
tions found in the spectral densities, the density of states 
and the conductivity, these effects are certainly not as 
spectacular. Mathematically this results from the fact 
that S z (k) involves an integration over energy w while 
it is this energy dependence which reflects most directly 
the details of the correlation effects. 

The in-plane spin texture is also changed, although the 
perpendicular spin momentum locking remains. In Fig. 8 
we show results for the momentum dependence (in units 
of k/kF) of S x -y(k) defined in Eq. (19). It gives the 
magnitude of the in-plane spin at momentum k. Without 
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interactions it would be equal to , „ „ „ „ t 

and is shown as the solid black curve for a gap of 20 meV 
and temperature T = 0.1K. It starts at zero for k = 0 
and rises monotonically with a sharp cut off at k/k F = 1 
at zero temperature. When the electron-phonon inter¬ 
action is included through Eq. (20) we get the dotted 
blue curve which follows closely the solid black curve 
with the important difference that it has tails beyond 
k = k F characteristic of correlation effects. The dashed 
red curve is for the same value of g but now the gap has 
been reduced to 10 meV. In this case the curve rises more 
sharply out of k = 0, flattens before showing a sharp drop 
at k = k F and the correlation tails beyond this drop off 
are now considerably larger. The correlation tails seen in 
S z (k) as well as in S x - y (k) are very similar to the much 
more familiar case22 of the momentum distribution in a 
fermion system with electron-phonon interaction. The 
probability of occupation of the state k at zero tempera¬ 
ture denoted by n(k) is given in Eq. (17) and in the bare 
band picture is a step function 0{k F — k). As shown in 
Fig. 9 long dashed black curve for T = 0.1K, g = 10 and 
a gap A = 20 meV, n(k) is considerably reduced from 
value one throughout the occupied states. It still has a 
finite discontinuous jump at k = k F with correlation tail 
for k > k F . The short red dashed curve is S z (k) previ¬ 
ously calculated and presented here for comparison with 
n(k). This shows the qualitatively similar effect of the 
electron-phonon interaction on these two quantities. 


V. CONCLUSIONS 


first function A(k,s,ui) can be measured directly in an¬ 
gular resolved photo-emission spectroscopy (ARPES). 
Its average over momentum k determines the electronic 
density of states N(u>) measured in scanning tunneling 
spectroscopy (STM). The second A z (k, s, oj) enters along 
with A(k, s, oj) the expressions for the dynamic longitu¬ 
dinal optical conductivity. It also determines the out 
of plane spin texture. Both A and A z are functions 
of oj and are encoded with sharp phonon structures at 
oj = ±c Je, ~g~ ue — A and — g — oje + A. These phonon 
structures manifest directly in the density of states and 
in a somewhat different way in the optical conductivity. 

The out of plane and in plane spin textures are modi¬ 
fied by the electron-phonon interaction. But these quan¬ 
tities involve an integration over energy of an overlap of 
spectral density and thermal factor. Consequently the 
manifestation of phonon-electron coupling in these quan¬ 
tities is more subtle and not as direct. Nevertheless im¬ 
portant corrections to a bare band picture arise. Just 
as the probability of occupation of a state of momentum 
k is reduced from one for k less than k F , and extended 
(in momentum) tails appear for k > k Fl the magnitude 
of S z (k) follows the same trend. A similar picture ap¬ 
plies for the magnitude of the in plane spin component. 
While the spin remains locked in the direction perpen¬ 
dicular to momentum, correlation tails appear beyond 
k = k F . While the calculations are for the explicit case 
of the electron-phonon interaction, they serve to illus¬ 
trate how many body interactions in general modify the 
bare band picture. 


We calculated the effect of electron-phonon coupling 
on a system of helical gapped Dirac fermions. A simple 
Holstein model with coupling to a single Einstein phonon 
of energy oje was used. For massive Dirac fermions two 
self energies need to be introduced. There is the fa¬ 
miliar quasiparticle renormalizarion which changes the 
bare band energies to dressed dispersion curves and pro¬ 
vides damping. In addition there is a second self en¬ 
ergy directly associated with modifications of the gap. 
It is energy (u>) dependent and complex with the real 
part directly modulating the magnitude of the bare band 
gap. Both these self energies enter the quasipaticle spec¬ 
tral density A(k,s,oj) with k momentum oj energy and 
s = ± giving conduction and valence band respectively. 
In direct analogy, a second spectral density A z (k,s,oj) 
associated with the gap channel is also introduced. The 
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